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We construct families of symmetric, antisymmetric, and asymmetric solitary modes in one- 
dimensional bichromatic lattices with the second-harmonic-generating (x ) nonlinearity concen- 
trated at a pair of sites placed at distance I. The lattice can be built as an array of optical waveguides. 
CN| ' Solutions are obtained in an implicit analytical form, which is made explicit in the case of adjacent 

nonlinear sites, I = 1. The stability is analyzed through the computation of eigenvalues for small 
■ perturbations, and verified by direct simulations. In the cascading limit, which corresponds to large 

1 mismatch q, the system becomes tantamount to the recently studied single-component lattice with 

two embedded sites carrying the cubic nonlinearity. The modes undergo qualitative changes with 
the variation of q. In particular, at I > 2, the symmetry-breaking bifurcation (SBB), which creates 
asymmetric states from symmetric ones, is supercritical and subcritical for small and large values 
0^ | of Q, respectively, while the bifurcation is always supercritical at Z = 1. In the experiment, the cor- 

responding change of the phase transition between the second and first kinds may be implemented 
by varying the mismatch, via the wavelength of the input beam. The existence threshold (minimum 
total power) for the symmetric modes vanishes exactly at q = 0, which suggests a possibility to 
create the solitary mode using low-power beams. The stability of solution families also changes with 
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I. INTRODUCTION AND THE MODEL 

ov 

O' 

■ The structure of bound states in linear systems follows the symmetry of the underlying potential, a commonly known 
example being wave functions of eigenstates in symmetric double- well potentials [l| . The addition of the self-attractive 
nonlinearity leads to a qualitative change of the situation, causing the transition from the symmetric ground state 
^vq \ to an asymmetric one, if the strength of the nonlinearity exceeds a critical value Q. This transition was studied in 
detail for Bose-Einstein condensates (BECs) loaded into double- well potentials Q . Experimentally, the transition was 
realized in BEC B and in nonlinear optics, where a double-well structure was induced in a photorefractive material 
!• 

The limit case of the double- well setting with a tall potential barrier between the wells corresponds to the dual-core 
system, such as optical fibers @ and Bragg gratings Q with the twin-core structure, and pairs of linearly coupled planar 
waveguides with the (second-harmonic-generating) intrinsic nonlinearity [||. The symmetry-breaking bifurcation 
(SBB), which destabilizes the symmetric ground state in nonlinear dual-core systems and gives rise to an asymmetric 
one, was first discovered in the discrete self-trapping model [9j. In nonlinear optics, the SBB was predicted for 
continuous- wave (spatially uniform) states [|| , and for solitons ip], Gil ; m the models of dual-core fibers and Bragg 
gratings. The SBB was also analyzed for matter-wave solitons in the BEC loaded into double-channel potential traps 

The self-focusing cubic nonlinearity (the Kerr term in optics) gives rise to the symmetry-breaking phase transition of 
solitons of the first kind (alias the subcritical SBB) in the twin-core system. In that case, the branches of asymmetric 
modes emerge as unstable ones, going backward and then stabilizing themselves at turning points, from which they 
continue the evolution in the forward direction [Tfij . The character of the phase transition alters to the second kind 
(i.e., the SBB type changes from sub- to supercritical) under the combined action of the self- focusing nonlinearity 
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and a periodic potential acting along the free axis (transverse to the double-well potential) fill llal . The supercritical 
SBB gives rise to stable asymmetric branches which immediately go in the forward direction [15j. The SBBs in the 
twin-core Bragg grating and x^ waveguides are of the forward type too @, H]- 

A noteworthy counterpart of the linear double-well potential is a pseudopotential [l6[ induced by a spatial mod- 
ulation of the nonlinearity coefficient in the form of a symmetric pair of sharp peaks. This configuration may be 
implemented in optics and BEC alike [13, EH . The ultimate form of the double-peak pseudopotential features the 
nonlinearity concentrated at two points, in the form of a symmetric pair of delta-functions (which may be approxi- 
mate d by narrow Gaussians) [19l - t21| . The SBB of solitons in the dual-core pseudopotentials belongs to the subcritical 
type dlUl. The SBB of localized modes was also recently investigated in the model of a linear waveguide with two 
narrow stripes embedded into it [22j (the solution for the mode pinned to a single stripe was found in Ref. [23]). 

For discrete solitons in dual-core lattices, the SBB was analyzed too, assuming the uniform transverse coupling 
between parallel chains (24|, or the coupling applied at a single pair of sites [2]|. The bifurcation is subcritical in the 
former case, and supercritical in the latter one. The simplest realizations of the SBB in discrete media are provided 
by linear lattices with a symmetric pair of embedded nonlinear sites (this setting was introduced in Refs. [26l. l2^|). 
or a symmetric pair of nonlinear elements side-coupled to the linear chain [28j]. In particular, the SBB for localized 



modes in the linear lattice with a pair of sites carrying the cubic nonlinearity was recently studied in Ref. [27| , where 
it was concluded that the bifurcation is, chiefly, of the subcritical type (except for the case when the nonlinear sites 
are separated by a single lattice spacing, see below). 

The objective of the present work is to investigate localized symmetric, asymmetric, and antisymmetric solitary 
modes in the bichromatic one-dimensional linear lattice with a pair of inserted x^-nonlinear sites. The lattice is 
described by the following equations: 

.dum + 1 ^ _|_ _ 2 Um ) + (S m0 + 5 m - i,o) u* m v m = 0, (1) 
dz 2 

H^P- + i (v m +x + u m _i - 2v m ) - qv m + i (8 m + 5 m _i fi )u 2 m = 0, (2) 
dz 2 2 

where u m and v m are complex amplitudes of the fundamental-frequency (FF) and second-harmonic (SH) fields at the 
m-th site, the constant of the linear coupling between adjacent sites is scaled to be 1, as well as the strength of the 
X^ nonlinearity inserted at sites m = and m = I, which are represented by the Kronecker's symbols <5 lrij o and 5 m —i o 
(i.e., integer I is the distance between the nonlinear sites), and q is the real mismatch parameter. The system can be 
implemented as an array of parallel optical waveguides, with z being the propagation distance along the waveguides, 
as previously done for a great variety of models supported by such quasi-discrete settings in nonlinear optics j29j . 
Arrayed x structures may be also built by means of the quasi-phase- matched technique [30j . In these contexts, 
discrete x^ systems were introduced in Ref. |3l| . 

To realize the present model, based on Eqs. (|T|), ([2"]). two selected cores in the waveguiding array can be made 
quadratically nonlinear by fabricating them of an appropriate material, or a pply ing the x^-inducing poling to this 
pai r [22| ] . In fact, the model with two Kerr- nonlinear sites, studied in Ref. [27|. corresponds to the cascading limit 
[32| of Eqs. ([T]), ([2]) for q — > +00. In this work, we demonstrate that the system with negative, zero, and positive 
values of the mismatch (q) opens new possibilities for the creation and manipulations of discrete solitary modes, in 
comparison with the cascading limit. In particular, we demonstrate that the character of the SBB can be switched 
from super- to subcritical by varying the mismatch, which also alters the stability of the modes and their existence 
thresholds. 

The paper is organized as follows. Analytical solutions for symmetric, asymmetric, and antisymmetric localized 
modes generated by Eqs. (UJ) , (J2]) in the infinite lattice are produced in Section 2. In the general case, the solution is 
implicit, while explicit solutions are obtained for the smallest separation between the nonlinear sites, 1 = 1. Numerical 
results, obtained for finite lattices, are reported in Section 3. The stability of the discrete modes is considered in that 
section too. The paper is concluded by Section 4. 

II. ANALYTICAL CONSIDERATIONS 
A. The general case 

Analytical solutions to Eqs. (TTJ), (|2|) for stationary modes with propagation constant k are sought for in the usual 
form, 

Akzjj 2ikz 



K%}= {e ikz U m ,e 2ik *V m }, (3) 
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which reduces Eqs. (HJ and ([2]) into the stationary equations for real discrete fields U m , V m 



kU„ 



- {Um+l + U m -1 — 2J7 TO ) + (<5 m ,0 + Sm-lfi) 



o, 



- (4fc + q) V m + - (Kn+l + V m -i - 2V m ) + - (S mfi + 5 m - h0 ) U 2 n = 0. 

At m < 0, an exact solution to Eqs. Q and ([5]) is obvious: 

U m = A\ exp (— Ki|m|) , V m = A 2 exp (-re 2 |m|) , 
where A\ and A2 are arbitrary amplitudes, and «i. 2 > are determined by the following relations: 

k = 2 sinh 2 («i/2) , 4fc + = 2 sinh 2 (« 2 /2) . 
Due to condition Ki t 2 > 0, the limit value of k which corresponds to «i = or k 2 — can be found from Eq. (|7]l : 



0, for g > 0, 
-g/4, for g < 0, 



the localized modes existing at k > kn m . Similarly, at m > I, the exact solution to Eqs. (01 and §5§ is 

U m = Ci exp (— Ki(m - 0) > Ki = C 2 exp (— « 2 (m - Z)) , 

and in the inner region, < m < I, the solution is constructed as 

U m = Bix exp (— Kim) + B 12 exp(— k% (Z — m)) , 
Kn = -B21 exp (— K 2 m) + i?22 exp (— « 2 (Z - m)) . 



(4) 
(5) 

(6) 
(7) 

(8) 

(9) 

(10) 



The conditions of the continuity of the discrete fields at points m = and m = I lead to linear relations between 
the amplitudes of the solutions in the outer and inner regions: 



Bn + B12 exp (-k-lI) = Ax, B n exp(—Kil) + Bi2 = C 1 , 
B21 + .B 22 exp(-/« 2 Z) = A 2 , B 2 \ exp (-k 2 Z) + B22 = C 2 . 

Finally, the nonlinear equations at sites m = and m = Z take the following form: 



(11) 



fc + 1 - i CXp (-Kl) 



-- [B n exp (-«!) + Su exp (- (/ - 1) Kl )] + A X A 2 



= 0, 



4k 



1 - 2 CX P(~ K2 ) 



Ao 



■- [B 2 i exp (-k 2 ) + S 22 exp (- (Z - 1) « a )] + l -A\ = 0, 



(12) 



fc + 1- -exp(-«i) 



Ci 



-- [B n exp (-«i (Z - 1)) + B u exp (-«i)] + CiC 2 



4fc + g + 1 - - exp (— Ka) 



C 2 



0. 



+ i [B 21 exp (-/ta (Z - 1)) + B 22 exp (-k 2 )] + ^C 2 = 0. 



(13) 



Amplitudes -611,12 and £?2i,22 can be eliminated in favor of A\^ and Ci i2 , making use of linear equations (|11[) , which 
leads to the system of four equations, (|12[) and (|13p . for four remaining unknown amplitudes, A\p and Ci j2 . Families 
of solutions for localized modes are characterized below by dependences of their powers (norms), 



m— +00 



N. 



{U,V} 



= V W 2 v 2 ) 



(14) 
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on the propagation constant, k, a dynamical invariant of Eqs. ([T]), ([2]) being the total power, 

N = Nu + AN V . (15) 

The solutions for the symmetric and antisymmetric modes are defined, respectively, by constraints A„ = C n , 
B n i = B n2 , or A n = (— l) n C n , B n \ — (—l) n B n 2, with n = 1,2. In the latter case, only the FF field is antisymmetric, 
while its SH counterpart remains symmetric; note also that the antisymmetric modes with even and odd I may be 
categorized, respectively, as ones of the on-site and inter-site types [231. Explicit solutions for the symmetric and 
antisymmetric modes can be obtained, in a simple approximate form, for < k, 4fc + q <C 1, when Eq. ([7]) yields 

Kt « V2k, k 2 « y/2(4k + q), (16) 

hence the modes are broad in this limit, according to Eq. ((6]). Further, the amplitudes and powers (|14j) of the 
symmetric and antisymmetric modes, found in the same limit, are 

symm: A\ w k 2 A 2 w (1/2) Ki k 2 , N v w k 2 /2, AT v w k\/ (4« 2 ) , 
antisymm: A 2 « k 2 A 2 « k 2 /Z, A^y w « 2 / (Z«i) , Ay « (Z 2 /t 2 )~ . (17) 

A corollary of this result is that, except for the case of q = 0, when k 2 2/ci, the total power of the modes never 
vanishes at ki j2 — > 0, hence there is a finite power threshold (a minimum value of the total power) necessary for the 
existence of the localized sta tes a t g ^ 0. In particular, for the symmetric modes with < q <C 1, the threshold is 
A r m i n = Nu{k = Ki = 0) w y/q/2, which, indeed, vanishes solely at = 0. 

On the other hand, in lattices of a finite size, the threshold power of the symmetric modes vanishes in the limit of 
Ki ~ v2k — > [i.e., when k 2 remains finite, provided that q is positive, see Eq. (|16[) ]. Indeed, in this limit case the 
top line in Eq. (|17[) implies the vanishing of both A\ and A 2 , while the nonvanishing of the threshold power in the 
respective infinite lattice is accounted for by the simultaneous divergence of the spatial width of the FF component, 
K7 1 — > oo. Obviously, the latter factor cannot compensate the vanishing of the amplitude in finite-size lattices. 



B. The case of I = 1: Explicit results 



The above analysis yields results in the implicit form, given by Eqs. ([Ti" l) -([l3 ]l . Explicit solut ions can be obtained 
for the smallest distance between the nonlinear sites, 1 = 1. In this case, one does not need to introduce solutions 
(|10[) for the inner layer, while the remaining equations for amplitudes A\_ 2 and C\_ 2 take the following form: 



k + 1- -exp(-Ki) 



4fc + q + 1 - - cxp (-K2) 



Ai + iCi + AiAa = 0, 



Ao 



2 C2+ 2^ 



0. 



(18) 
(19) 



k + 1 - - CXp (-Kl) 



Ak + q + 1 - - exp (-k 2 ) 



d + -Ai + C X C 2 

C 2 + \a 2 + \c\ 



0, 
0. 



(20) 
(21) 



The symmetric solution, with A\ — C\ and A 2 = C 2 , can be easily found from here: 



, 11 , \ 
A 2 = k+ - - -exp(-/ci). 



4fc + (?+ - - icxp(-K 2 ) 



.4, 



(22) 



The search for the SBB, i.e., solutions with infinitely small A\ — C\ and A 2 — C 2 , yields an equation for the bifurcation 
point: A\ = 4fc + q + 1 — | cxp (— k 2 ). Substituting A 2 from Eq. (|22|) . it takes the following form: 



4fc + g+ i - iexp(-« 2 ) 



[2fc — cxp (—Ki)] = 1. 



(23) 
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Further analysis demonstrates that, at I = 1, the SBB is supercritical at all values of q, including the cascading 
limit, q — > oo. The latter fact implies that the SBB is also supercritical in the model with the cubic nonlinearity at 
two adjacent sites. Indeed, this can be checked to be correct in the cubic model, while for alH > 2 the corresponding 



SBB is subcritical [27[ (the character of the SBB for I = 1, super- or subcritical, was not considered in Ref. [271 ]) 
It is also possible to find antisymmetric solutions, with A\ = — Ci, = C2: 



3 1 

2-2 eXp( " 



4k 



1 



CXp (— K-i 



A 2 



(24) 



A point where the antisymmetry could be spontaneously broken, similar to the SBB of the symmetric modes, is 
determined by the condition that solutions with infinitely small A\ + C\ and A2 — C2 emerge. After a simple analysis, 
this condition leads to the following equation: 



4k 



1 1 < 



K2) 



[2k + 2 - exp(-Ki) 



(25) 



cf. Eq. (|23[) . Replacing here 4k + q and k in the square-bracket combinations by their expressions in terms of K2 
and K\ given by Eqs. ([Tj) , it is easy to check that Eq. ([25]) can never be satisfied, unlike Eq. (|23| (the left-hand side 
always takes values > 1), hence the antisymmetric modes do not undergo the bifurcation. 



III. NUMERICAL RESULTS 
A. Zero mismatch (q — 0) 

First, we present families of numerical solutions of Eqs. (|4|) and ([5]) with a finite number of sites, Af = 91, obtained 
for q = and I = 5. In Fig. [TJ the families of symmetric, asymmetric, and antisymmetric modes are represented 
by the Njjy{k) curves for the FF and SH components. These curves completely overlap with those predicted by the 
implicit analytical solutions, in the form of Eqs. (p~T|) - (fT3|) . 

In Fig. Ufa), the family of symmetric solutions features the SBB at a finite value of fc, and an asymptotic behavior 
with Njjy — > at k — > (i.e., the vanishing of the threshold power), in accordance with the analytical approximation 
based on Eqs. (fll)|) and (fTT|) . This result is of obvious interest to the experiment, suggesting the possibility to create 
the nonlinear guided modes, using low-power input beams. On the other hand, Fig. |T](b) demonstrates a completely 
different behavior of the norms of the antisymmetric solutions at k — > 0: While Nu approaches a finite value in this 
limit, Ny diverges at k — > 0, also in accordance with Eqs. (|T7|) and (|16[) (in terms of the numerical results, both the 
divergence of Ny and maintaining the finite value of Njj arc limited by the finite size of the lattice). 

The stability of the localized modes of the different types, which is also indicated in Fig. [T] was established by 
solving the linear eigenvalue problem for small perturbations, and verified through direct simulations of perturbed 
evolution in the framework of Eqs. (HJ) , (BJ) - 

Typical profiles of modes featuring the different types of the symmetry are shown in Fig. ® an d the (in) stability of 
these modes is illustrated in the first three columns of Figs. [3]by means of direct simulations. As might be expected, 
the instability of the symmetric localized modes past the SBB point (fcgBB ~ 0.015) transforms it into an excited state 
close to a stable asymmetric mode. The simulations also corroborate the change of the stability of the antisymmetric 
modes [at k 0.075 in Fig. [T^b)], which is predicted by the computation of eigenvalues for small perturbations. In 
particular, the fourth column of Fig. [3] demonstrates an example of stable antisymmetric modes. 

To highlight the character of the SBB, we use natural measures of the asymmetry of the solutions taken in the form 
of Eqs. (j5])- (fTU)) . namely. 



A 2 - C 2 A 2 - r 2 

Al + CV A 2 + C? 



Vu^-^z, Vv^-rt-rtl- (26) 



For g = 0, the asymmetries are plotted versus k and the total power (norm), defined as per Eq. (fi"5j) . in Fig. [H which 
demonstrates that the bifurcation is supercritical in this case. 



B. Negative mismatch (q < 0) 



A behavior different from that at q = is observed at q < 0. As seen in Fig. not only the antisymmetric 
solutions but also symmetric ones feature the divergence of the norm of the l/-field as k approaches the corresponding 
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FIG. 1: (Color online) The powers (norms) of the FF (U) and SH (V) components of one-dimensional localized modes versus 
the propagation constant, as produced by the implicit analytical solution for the infinite lattice, based on Eqs. (|lip - (|13| l. and 
by the numerical solution for the finite lattice with two nonlinear sites, separated by distance I — 5, with zero mismatch, q — 0. 
The analytical and numerical curves completely coincide. Families of symmetric and asymmetric solutions, and antisymmetric 
ones, are shown, severally, in panels (a) and (b). Solid and dashed lines correspond to the stable and unstable branches, 
respectively. 




FIG. 2: (Color online) Examples of analytical (solid line) and numerical (points) solutions for symmetric (unstable, in this 
case), asymmetric (stable), and antisymmetric (unstable) modes, found for k = 0.02, / = 5 and q = 0. Black and red colors 
show the U- and V^-fields, respectively. 

limit value fcn m = —q/^, see Eq. (j8|). This feature can also be easily explained by the above analysis: In the limit of 
k + q/4 —> 0, Eq. ([7]) yields K2 ~ \J1 (4fc + q), hence the respective width ~ k^ 1 diverges, while K\ remains finite [cf. 
Eq. (|16j) ]. Further, Eqs. (|li p -(|13 j) yield, in this limit, finite A2 for both symmetric and antisymmetric branches, while 
the respective value of A\ is vanishing as per A\ ps K2A2 cf. Eq. (fT7|) ]. hence the corresponding norm Njj vanishes 
too, while Ny pa n^ 1 A\ diverges. 

The linear-stability analysis demonstrates that the symmetric branch is stable up to the SBB point (/csbb ~ 0.0514). 
Past the bifurcation point, the asymmetric branch appears in an unstable form, and with the further increase of k 
it becomes stable at k ~ 0.064. Direct numerical simulations displayed in Fig. [5] confirm the predictions of the 
linear-stability analysis. In particular, the unstable asymmetric and antisymmetric modes are transformed by the 
perturbed evolution into breathers. 

Finally, the SBB in the case of q < is of the supercritical type, similar to that displayed in Fig. 0]for q = 0. 

C. Positive mismatch (q > 0) 

The increase of mismatch q to large positive values leads to a transition between the supercritical and subcritical 
types of the SBB. In particular, Fig. [7] shows that the Njj(k) branch of the asymmetric solutions passes a shallow 
minimum right after the bifurcation point, which is a signature of a subcritical bifurcation, while the bifurcation 
in Fig. Ufa) was supercritical. The change of the type of the SBB at q > q n from super- to subcritical is further 
illustrated by Fig. [5] It is worth noting that the bifurcation is subcritical only in terms of the FF field, U. 

For all I > 2 it is possible to find a boundary value qo, such that the SBB is of the super- and subcritical types at 
q < qo and q > qo, respectively Numerical analysis of Eqs. (|12p and (|13p demonstrates that go rapidly grows with the 
decrease of I: qo(l = 5) ps 0.27, qo(l = 3) pa 1.86, and qo(l = 2) ps 4.9. As said above, the SBB keeps its supercritical 
character at all values of q for I = 1. The possibility to switch between the different types of the SBB, i.e., between 
the phase transitions of the second and first kinds by means of the mismatch, may be readily implemented in the 
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FIG. 3: (Color online) The first three columns display density plots representing the perturbed evolution of symmetric (unsta- 
ble), asymmetric (stable) and antisymmetric (unstable) solutions from Fig. [2] The fourth column displays the density plot of 
the stable antisymmetric solution at k = 0.08. 




FIG. 4: (Color online) The dependence of asymmetry measures Q26p on propagation constant k and total power N in the case 
of the supercritical bifurcation, at I = 5 and q = 0. Black and red colors show the U- and V-fields, respectively. 



experiment, as the mismatch may be varied by means of the wavelength of the input beam. 

The trend of the transition to the subcritical bifurcation with the increase of q complies with the fact that, as 
mentioned above, at large q the cascading approximation applies to Eqs. (|4]), ([5]), making them asymptotically 
tantamount to the equation for the single-component (FF) lattice with the two sites carrying; the cubic self-focusing 
nonlinearity. In the latter case, the SBB for the localized modes is subcritical for all / > 2 [271 ]. 

The results of the analysis performed at all values of mismatch q, negative, zero, and positive, are summarized in 
Fig. [9l where the critical value of the propagation constant at which the SBB happens, fcsBB, is plotted as a function 
of q at different fixed values of distance I between the nonlinear sites. Note that, at large I, the bifurcation point 
&sbb is very close to the limit value ku m given by Eq. ([8]), which is simply explained by the fact that large I implies 
an exponentially weak coupling between the two nonlinear sites, hence the stability margin of the symmetric state is 
vanishingly small. On the other hand, at large q > the curves approach constant values, which correspond, in the 
cascading limit, to the system with the cubic nonlinearity. 



IV. CONCLUSION 



The objective of this work is to develop the analysis of the stability and spontaneous symmetry breaking of discrete 
solitary modes in one-dimensional lattices with the quadratic nonlinearity, in the most basic setting with the nonlin- 
earity applied at a two sites embedded into the linear host lattice, with distance I between them. The system can be 
readily implemented as an optical waveguiding array. The solutions for symmetric, asymmetric and antisymmetric 
modes are obtained in the implicit analytical form, which becomes explicit for 1 = 1. The stability of the modes was 
investigated through the computation of eigenvalues for small perturbations, and checked against direct simulations 
of the evolution of perturbed modes. The analysis has demonstrated that properties of the modes undergo essen- 
tial evolution with the change of the mismatch, q: The SBB (symmetry-breaking bifurcation) changes from super- 
into subcritical (except for the case of I = 1), the existence threshold for the symmetric localized modes vanishes 
exactly at the point of q = 0, and the stability of the branches changes too. These results suggest a straightforward 
implementation in the experiment. 
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FIG. 5: (Color online) Panels (a) and (b) show the same as in Fig. [TJ but for q = —0.2. Panel (c) is a blowup of the bifurcation 
region for the (7-field. 
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FIG. 6: (Color online) The evolution of symmetric, asymmetric and antisymmetric unstable modes at k = 0.054, for I — 5 and 
q = -0.2. 



The analysis reported in this work can be naturally extended in other directions. In particular, following Refs. [27[ 
and [2l| . it may be interesting to consider a finite ring-shaped lattice with the nonlinear sites set at diametrically 
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FIG. 7: (Color online) The SBB diagram of the localized modes, as produced by the analytical solution for the infinite lattice 
with two nonlinear sites, separated by distance I — 5, and mismatch q — 5 (the numerical solution for the same case is 
indistinguishable from the analytical one). Solid and dashed lines correspond to stable and unstable branches, respectively. 




FIG. 8: (Color online) The same as in Fig. SI but for q = 5, with the respective bifurcation being subcritical. 



opposite points. A challenging problem is to generalize the analysis for two-dimensional lattices. 



[1] D. Landau and E. M. Lifshitz, Quantum Mechanics (Moscow: Nauka Publishers, 1974). 

[2] E. A. Ostrovskaya, Y. S. Kivshar, M. Lisak, B. Hall, F. Cattani, and D. Anderson, Phys. Rev. A 61, 031601 (2000); R. 

D'Agosta, B. A. Malomed, C. Presilla, Phys. Lett. A 275, 424 (2000); R. K. Jackson and M. I. Weinstein, J. Stat. Phys. 

116, 881 (2004); D. Ananikian and T. Bergeman, Phys. Rev. A 73, 013604 (2006); E. W. Kirr, P. G. Kevrekidis, E. 

Shlizerman, and M. I. Weinstein, SIAM J. Math. Anal. 40, 566 (2008). 
[3] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, Phys. Rev. A 55, 4318 (1997); A. Smerzi, S. Fantoni, S. Giovanazzi, 

and S. R. Shenoy, Phys. Rev. Lett. 79, 4950 (1997); S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, Phys. Rev. 

A 59, 620 (1999); K. W. Mahmud, H. Perry, and W. P. Reinhardt, Phys. Rev. A 71, 023615 (2005); E. Infeld, P. Zih, 

J. Gocalek, and M. Trippenbach, Phys. Rev. E 74, 026610 (2006); G. Theocharis, P. G. Kevrekidis, D. J. Frantzeskakis, 

and P. Schmelcher, Phys. Rev. E 74, 056608 (2006); G. L. Alfimov and D. A. Zezyulin, Nonlinearity 20, 2075 (2007); C. 

Wang, P. G. Kevrekidis, N. Whitaker, and B. A. Malomed, Physica D 237, 2922 (2008). 
[4] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M. K. Oberthaler, Phys. Rev. Lett. 95, 010402 (2005); R. 

Gati, M. Albiez, J. Foiling, B. Hemmerling, and M. K. Oberthaler, Appl. Phys. B 82, 207 (2006). 
[5] P. G. Kevrekidis, Z. Chen, B. A. Malomed, D. J. Frantzeskakis, and M. I. Weinstein, Phys. Lett. A 340, 275 (2005). 
[6] A. W. Snyder, D. J. Mitchell, L. Poladian, D. R. Rowland, and Y. Chen, J. Opt. Soc. Am. B 8, 2101 (1991). 
[7] W. C. K. Mak, B. A. Malomed, and P. L. Chu, J. Opt. Soc. Am. B 15, 1685 (1998). 
[8] W. C. K. Mak, B. A. Malomed, and P. L. Chu, Phys. Rev. E 55, 6134 (1997); ibid. 57, 1092 (1998). 
[9] J. C. Eilbeck, P. S. Lomdahl, and A. C. Scott, Physica D 16, 318 (1985). 
[10] C. Pare and M. Florjariczyk, Phys. Rev. A 41, 6287 (1990); A. I. Maimistov, Kvant. Elektron. 18, 758 [Sov. J. Quantum 

Electron. 21, 687 (1991)]; N. Akhmediev and A. Ankiewicz, Phys. Rev. Lett. 70, 2395 (1993); P. L. Chu, B. A. Malomed, 

and G. D. Peng, J. Opt. Soc. A B 10, 1379 (1993); B. A. Malomed, in: Progr. Optics 43, 71 (E. Wolf, editor: North 

Holland, Amsterdam, 2002). 



10 




FIG. 9: (Color online) The critical value of the propagation constant, £;sbb, at which the symmetry-breaking bifurcation occurs, 
as a function of mismatch q, at fixed I = 1, 2 and 5. Open squares for / = 1 correspond to analytical prediction (|23|l . The 
dashed line corresponds to fc m , as determined by Eq. ([8]). 



[11] A. Gubeskys and B. A. Malomed, Phys. Rev. A 75, 063602 (2007). 

[12] M. Matuszewski, B. A. Malomed, and M. Trippenbach, Phys. Rev. A 75, 063621 (2007). 

[13] M. Trippenbach, E. Infeld, J. Gocalek, M. Matuszewski, M. Oberthaler, and B. A. Malomed, Phys. Rev. A 78, 013603 

(2008); N. V. Hung, M. Trippenbach, and B. A. Malomed, ibid. Phys. Rev. A 84, 053618 (2011). 
[14] L. Salasnich, B. A. Malomed, and F. Toigo, Phys. Rev. A 81, 045603 (2010). 

[15] G. Iooss and D. D. Joseph, Elementary Stability Bifurcation Theory (Springer- Verlag: New York, 1980). 

[16] W. A. Harrison, Pseudopotentials in the Theory of Metals (Benjamin: New York, 1966). 

[17] Y. V. Kartashov, B. A. Malomed, and L. Torner, Rev. Mod. Phys. 83, 247 (2011). 

[18] L. C. Qian, M. L. Wall, S. Zhang, Z. Zhou, and H. Pu, Phys. Rev. A 77, 013611 (2008). 

[19] T. Mayteevarunyoo, B. A. Malomed, and G. Dong, Phys. Rev. A 78, 053601 (2008). 

[20] N. Dror and B. A. Malomed, Phys. Rev. A 83, 033828 (2011). 

[21] X.-F. Zhou, S.-L. Zhang, Z.-W. Zhou, B. A. Malomed, and H. Pu, Phys. Rev. A 85, 023603 (2012). 
[22] A. Shapira, N. Voloch-Bloch, B. A. Malomed, and A. Arie, J. Opt. Soc. Am. B 28, 148 (2011). 
[23] A. A. Sukhorukov, Y. S. Kivshar, and O. Bang, Phys. Rev. E 60, R41 (1999). 

[24] G. Herring, P. G. Kevrekidis, B. A. Malomed, R. Carretero-Gonzalez, and D. J. Frantzeskakis, Phys. Rev. E 76, 066606 
(2007). 

[25] Lj. Hadzievski, G. Gligoric, A. Maluckov, and B. A. Malomed, Phys. Rev. A 82, 033806 (2010). 

[26] M. I. Molina and G. P. Tsironis, Phys. Rev. B 47,15330 (1993); B. C. Gupta and K. Kundu, Phys. Rev. B 55, 894 (1997). 
[27] V. A. Brazhnyi and B.A. Malomed, Phys. Rev. A 83, 053844 (2011). 

[28] B. Maes, M. Soljacic, J. D. Joannopoulos, P. Bienstman, R. Baets, S.-P. Gorza and M. Haelterman, Opt. Express, 14, 
10678 (2006); E. N. Bulgakov and A. F. Sadreev, Phys. Rev. B 81, 115128 (2010); E. Bulgakov, A. Sadreev, and K. N. 
Pichugin, ibid. B 83, 045109 (2011). 

[29] F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. Assanto, M. Segev, and Y. Silberberg, Phys. Rep. 463, 1 (2008). 

[30] C. B. Clausen, P. L. Christiansen, L. Torner, and Y. B. Gaididei, Phys. Rev. E 60, R5064 (1999). 

[31] O. Bang, P. L. Christiansen, and C. B. Clausen, Phys. Rev. E 56, 7257 (1997); T. Peschel. U. Peschel and F. Lederer, ibid. 

E 57, 1127 (1998); S. Darmanyan, A. Kobyakov, and F. Lederer, ibid. E 57, 2344 (1998); S. Darmanyan, A. Kamchatnov, 

and F. Lederer, ibid. E 58, R4120 (1998). 
[32] G. I. Stegeman, G. J. Hagan, and L. Torner, Opt. Quantum Electron. 28, 1691 (1996). 



